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Abstract 

We show that an excehent approximation to the exact quantum solution of the ground state of 
the Tavis-Cummings model is obtained by means of a semi-classical projected state. This state has 
an analytical form in terms of the model parameters and, in contrast to the exact quantum state, 
it allows for an analytical calculation of the expectation values of field and matter observables, 
entanglement entropy between field and matter, squeezing parameter, and population probability 
distributions. The fidelity between this projected state and the exact quantum ground state is 
very close to 1, except for the region of classical phase transitions. We compare the analytical 
results with those of the exact solution obtained through the direct Hamiltonian diagonalization 
as a function of the atomic separation energy and the matter-field coupling. 

PACS numbers: 42.50.Ct, 03.65.Fd, 64.70.Tg 
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I. INTRODUCTION 



The progress in the technology of trapped-atom lasers and the cavity QED experiments 
has resumed the interest of researchers in the basic models that describe the interaction 
between quantized radiation and atoms. The Dicke Model [l] (DM) describes the inter- 
action of a quantized radiation field with a sample of two-level atoms located within a 
distance smaller than the wavelength of the radiation. Dicke realized that under certain con- 
ditions a gas of radiating molecules shows a collective behavior called superradiance. This 
phenomenon was observed experimentally in optically pumped HF gas [2|]. The simplest 
case = 1 in the rotating-wave approximation is known as the Jaynes-Cummings model 
(JCM) j^. Many theoretical predictions of this latter model, such as the existence of col- 



lapse and revivals in the Rabi oscillations 



J], the formation of macroscopic quantum states, 



or measures of entanglement associated with spin-squeezed states, have been confirmed, and 
many experimental studies of Rydberg atoms with very large principal quantum number 
within single-mode cavities have been observed [5!]. It is well known that the JCM has 
served as a guide to understand several quantum optics phenomena. While this considers 
the rotating wave approximation in order to discard the non-conserving energy terms, the 
model is widely applicable: in 1991 it was proposed [61] that by using circularly polarized 
light together with atomic selection rules the non-conserving energy terms can be elimi- 
nated; additionally, the diamagnetic term can be taken into account by properly changing 
the frequency of the field mode as a function of the coupling interaction strength: the 
renormalization of the field mode frequency can be understood by making a Bogoliubov 
transformation on the field part of the Dicke Hamiltonian. This implies that, at least in 
some cases, the results of the JCM can be valid beyond the rotating wave approximation. 

The JCM, which might well be called the standard model of quantum optics, has been 
studied intensively for more than forty years. An experimental example is provided by Aoki 
et al. 71]. When there are many atoms interacting with a single mode quantized radiation 
field of one and the same cavity, the exact solution can be obtained by means of the so- 
called Tavis-Cummings model |8| (TCM), and its recent generalizations 9[]. This model 
predicts the collective A^-atom interaction strength to be 7Ar = y/N gi, where gi is the dipole 
coupling strength of each individual atom i. The TCM also has been considered to describe 
cavity QED with a Bose-Einstein condensate 10| . Exact solutions of the generalized Tavis- 
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Cummings model have been found by means of quantum inverse methods, an algebraic 



procedure based on finding sets of solutions to the Bethe equations 



Other algebraic 



methods solve the eigenvalue problem of the TCM using polynomially deformed algebras [1 
where analytical expressions may be found up to third order in a specific perturbation theory. 

Since the presence of superradiant phase transitions in the DM were established [3| , there 
have been several contributions simplifying the original computation and others that have 



found also the presence of phase transitions in generalized Dicke models 



]J,ll5|,ll6i. In par- 



ticular, one of the contributions showed the existence of phase transitions in a Dicke Hamil- 
tonian which includes the counter-rotating atom-field interaction, although the necessary 
condition between the coupling parameter and the field mode frequency was modified [l^. 

In 1975 it was shown that the superradiant phase transition was due to the absence of the 
diamagnetic term in the Hamiltonian describing the N two-level atoms interacting with a one 
mode electromagnetic field l3] • It established, through the Thomas- Reiche-Kuhn sum rule. 



that the phase transition cannot be reached because it wou 
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d place contradictory bounds 



2m\ have established that gauge 



to the parameters of the model. Other authors 
invariance requires the presence of the diamagnetic and counter-rotating terms; moreover, 
the sum rule can be derived from charge-current conservation. For this reason, a criterion 
was seeked to establish the validity of Dicke models and determine whether a superradiant 
phase transition can occur. It was shown that two-level atoms which make electric dipole 
transitions cannot exhibit superradiant phase transitions, although their result does not rule 
out phase transitions based on magnetic dipole interactions. In the last decade, the work 
of Crisp was extended to many two-level atoms and it was concluded that in the strong 
coupling regime, in the thermodynamical limit, it is not possible to have superradiant phase 



transitions 



2l|. 



It is however known that, for a one electron atom interacting with a non-quantized 
electromagnetic field (that is, in the semiclassical approximation), a unitary transformation 
of the form f/(r, t) = exp [j^r ■ A(0, t) j can be performed to eliminate the terms depending 
on the electromagnetic vector potential, giving rise to an electric dipole-field interaction. 
The same procedure can be applied to the quantized electromagnetic field, the difference 
being that we are dealing with operators and there is an extra term, due to the commutation 



relations, which represents a dipole-dipole interaction 



22l |. In 2001, similar arguments on 



the interaction between radiation and molecular dipole moments were given [23[, where 
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they use a gauge invariant formulation of the molecular electric dipole-photon interactions 
which allows for superradiant phase transitions as in the Hepp-Lieb formulation. They 
show that the dipole-field interaction is strictly linear in the electric field, with an extra 
self-interaction term. Thus, there are no quadratic terms, and therefore a physical system 
displaying superradiant phase transitions should exist. The difference with previous gauge 
invariant formulations is due to the time dependence of the unitary transformation. 

The considerations above render likely the physical reality of systems presenting superra- 
diant phase transitions, and being described by a Hamiltonian of the type used to describe 



linear dipole-photon interactions. In particular, there are physica 



systems where interac- 



tions between bosonic degrees of freedom with a collective spin [24, l25|] are relevant. Then, 
the recently presented method |2Q] which allows a simple and elegant determination of the 
stability properties and phase transitions for a finite number of particles, is of importance. 

In this work we obtain an excellent approximation to the superradiant states of the ground 
state of the Tavis-Cummings Hamiltonian, which admits analytical expressions for both, field 
and matter observables, including the entanglement entropy between field and matter, the 
squeezing parameter, and the population probability distribution. The fitness between this 
approximation and the exact quantum solution is measured through the evaluation of the 
fidelity parametercitewooters. 

Section II below establishes the Tavis-Cummings Hamiltonian in terms of the constant of 
motion, and shows that the variational tensorial product of coherent states has difficulties 
in describing, for example, the photon number fluctuations. It is argued that restoring the 
symmetry by projecting the variational state to a given value of the constant of motion, a 
much better (and analytical) approximation is obtained. This projected state is justified 
and built. In Section HI, by means of the overlap of the projected state, the calculation of 
the expectation values for field and matter observables is done analytically for an arbitrary 
value of the constant of motion. The entanglement entropy between the matter and field, 
the squeezing parameter of the state and the probability distributions of photons and atoms 
are also calculated. In Section IV the selection of the appropriate value of the constant of 
motion A is discussed. Section V compares the results of Section HI with the corresponding 
ones for the exact quantum state, and Section VI draws some conclusions. 



4 



II. PROJECTED STATE 



The Tavis-Cummings model Hamiltonian for identical 2- level systems (e.g., atoms) 
immersed in an electromagnetic field, is given by 



Htcm = ^^f^'^^ + ^^aJz + -^{^'^ J- + ^J+) , (1) 



where up is the field frequency, uja the atomic energy-level difference, and the dipole cou- 
pling. The operators a, a^, denote the one- mode annihilation and creation photon operators, 
Jz the atomic relative population operator, and J± the atomic transition operators. 
It is immediate that this Hamiltonian commutes with the operator 

A = a% + . (2) 

It is then convenient to rewrite it by introducing a detunning parameter A = — u^, and 
by dividing it hj up (which can be thought of as the natural unit of frequency) and by the 
total number of particles, having in this way an intensive Hamiltonian operator 

where A = 1 - = 1 _ and 7 = ^. 



Following the variational procedure indicated in 2o|, it was found that the direct product 

n n 

of coherent states: Heisenberg-Weyl for the photon part \a) [51, |28| and SU{2) or spin for 
the matter \() [29|, i.e., \a, () = \a) ® IC), is a good approximation to the exact quantum 
solution of the Tavis-Cummings Hamiltonian, in spite of the fact that the overlap of the two 
states is small. This approximation describes very well the expectation values of the atomic 
observables and most of the field observables. However, particular difficulties are found in 
the description of the photon number fluctuations (An)^, as shown in the left of Fig.([T]), 
and for several observables when the detuning parameter A takes values greater than 1. 

The first of these problems arises because the variational state has contributions from all 
the eigenvalues A of the operator A, violating its conservation symmetry. It is then natural 
to propose as ground state a state with the symmetry restored by projecting the variational 
state to a given value of A, namely, the one that minimizes the classical energy obtained from 
the variational procedure and approximated to the closest integer or half-integer according 
to whether is even or odd, respectively. The exact procedure is presented in Section IV 
below. 
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FIG. 1: Squared fluctuation of the photon number operator n, for the ground state, as a function 
of the interaction strength 7, for a detuning parameter A = 0.2 and = 20 atoms. At left, 
the continuous line shows the coherent state approximation to the exact solution (discrete line). 
The differences are large even when normalized by the number of atoms. At right the lower curve 
shows the approximation we can obtain using a projected state, as explained in this paper. Note 
the different scales for both plots. 



By projecting the state to one (as yet arbitrary) value of A we obtain 



1^) 



|0) ® \j, -3) : 
'i/=max[0, A— j] 

|0) ® |j, 3) , 



> 7 ; 
\^a\ < 7^ ; 
< -7^ ■ 



(4) 



In this expression we have used j = -j, and for \uja\ < 7 the state is not normalized, with 



c 



A^7 



UJa 

T 



(5) 



It has the advantage that it only contains one eigenvalue of the constant of motion, it is an 
analytic solution, reproduces well the matter and field observables (cf. Fig ([T]) at right and 
note the scale difference with the figure at left), and the overlap with the exact quantum 
solution is very close to one {vide infra). 

The conditions on the parameters uja and 7 in Eq.(jl]) are associated to the existence of a 



separatrix in the model which defines the phase transitions in the semiclassical solution 26|. 

To build the projected state Eq.(jl]) we use the semiclassical procedure, which consists 
of the calculation of the expectation value of the Hamiltonian with respect to the tensorial 
product of coherent states \a) ® \x), the determination of the minima, and the use of the 



catastrophe formalism 



30| to find the stability properties |26|. We find that the critical 



points for a minimum are given by 



qc = -vi7 sin^'c cos(^c , 
Pc = sin 9c sin 0^ , 

where we have written a = + ip) in terms of the expectation values of the quadratures 
of the field, and x = e*'^tan(|), with {9,(j)) denoting a point in the unit Bloch sphere; at 
their critical values 6c, (pcy we have 26 1 



mmima: < 



0c = O, 

ec = 7T. 



Or 



arccos 



Eq 



N UJA 
2 ' 

2 ' 

, Af. 



Ac = -j 
Ac = 3 

j 



for uja > 7^ 
for uja < —7^ (6) 
for \uja\ < 7^ 



^^13 — V 4^2 ) '>C J 

This expression shows the minima critical points, the energy Eq, the constant of motion 
Ac, and the conditions in the parameter space to guarantee that they constitute an energy 
minimum. The energy surface of the Tavis-Cummings model is 0-unstable, for which reason 
0c can be taken arbitrarily. 

The expression for the trial state |a) |x) that minimizes the energy surface takes the 
following form 26l |: 



North Pole {uja > 1^): \^np) = |0) ® |j, -3) 
South Pole {ua < -7^): l^.p) = |0) ® |j, j) 



hj +00 



Parallels {\uja\ < 7^): li'par) = An,M) ® |j, m) 

m=—j u=0 

where we have defined the expansion coefficients 



(7) 
(8) 

(9) 



2j 
j + m 



1/2 



exp 



3T 
4 



UJ 



2 272 



1 - J + z (j + m - z/) 



1 



lua 



' i/\ \z ^7"/ v2 27^ 

One usually would perform this sum by selecting a value of u, making the sum over m, and 
then proceeding with the following value of z/, etc., until we reach some type of convergence. 
In the TCM model X = m + u is a conserved quantity. By replacing m by A — z/, we can 
write 26| 

+00 X+j 
A=— J j/=max(0,A— j) 



(10) 



The eigenstates of the North Pole ([7]), South Pole ([8]), and Parallel region (fTOl) . for a given 
A, justify the projected state established in (jl]) with the unnormalized state for \uja\ < 7^ 
defined by 

where we have retained only terms depending on the number of photons v from the coefficient 



III. EXPECTATION VALUES OF FIELD AND MATTER OBSERVABLES 

The effect of the transformation generated by the constant of motion exp[z(y9A] in the 
plane of quadratures of the electromagnetic field is a counterclockwise rotation by an angle 

while for the matter observables Jx and Jy it is a clockwise rotation by an angle ip along 
the 2;-axis. This transformation then leaves invariant each term in the TCM Hamiltonian. 
This gives our trial states their 0-invariance. Calculating for any value of allows one to 
recover the result for any other value of (j) through the rotation above. Without loss of 
generality, then, we choose = in what follows. 

To determine the expectation values of the observables of the system, in the Parallels 
region, it is useful to calculate the overlap 



A+i 



2j \ (C'C)^ 



;/=max[0, ^ ' 

Considering this expression as a truncated expansion of the Hypergeometric Confiuent func- 
tion we obtain 

{C'ICl = I • (13) 

where i^"(x) denotes tha associated Laguerre polynomials, which are defined as [3ll] 

Ll{x) = ^-^ ^ (e-^ x"+") , n > , 
^ fix" ^ ^ ~ 

L'-ni^) = nL°(x)-(n + a)L^,(x) , n>l. 

From expression (jlj) for the projected state, it is immediate that the probability of finding 
z/ photons depends of the relative values of uja and 7, i.e., 

8 



for \ua\ < 7^, while Vi, = Si^^q outside that region. The probabihty of finding Ue excited 
atoms can be obtained from the previous expression by replacing u X + j — rie- 

Defining i] = (^^, and the overlap of the projected state hjY = {CIC}? we can calculate the 
expectation value of the photon number operator n = a^a and its corresponding fluctuations 
squared (An)^ : 

(n) = J2 yV, = r]—\YiY , (15) 

i/=max[0, A— j] 
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(An)^ = [r^Yr,) " ^^^^ 

Note that for A = —j the overlap is y = 1. On the other hand, for X = j the overlap 
is y = L2j(— C^), which only equals 1 when uja = —7^- At the North and South Poles, 
therefore, the expectation values for n and (Ari)^ can be obtained by assuming Y = 1 there 
and the expressions (ITS!) . (fT6l) to be valid for all values of the parameters 7 and uja- Analytic 
expressions for (n) and (Ari)^ in the Parallels region are given in the Appendix. 

The expectation values of other matter and field observables can be determined in terms 
of the equations above. 

The quadrature components of the electromagnetic field are given by 

^ = ^ (a^ + a) , p = (at - a) ; 

as they change the value of A for any state, their expectation values with respect to the 
projected state are equal to zero. Their corresponding fluctuations are therefore 

if) = iP') = {n) + l. (17) 

For the expectation values of matter observables we have (J^.) = (Jy) = for the same 
reason as above, and others take the form 

(J,) = X-{n) , (18) 
{AJ^f = (Anf , (19) 
(J!) = {Jy) = lj{j + l)-l{^-{n)f (20) 

Finally, the expectation values for the transition operators a^J. and a J+ are given by 

(atj_) = (aJ+)=C(j + A-(n)) . (21) 
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IV. DETERMINATION OF THE CONSTANT OF THE MOTION 



The expectation value of the Hamiltonian with respect to the projected state for an 
arbitrary value of A is 



^''^ " icier 



+ A)c + ^ 



A 



27 



C 



(n) (22) 



When A = — j we have (n) = and this expression simplifies to (H) = — | (1 — A). By 
substituting the expression of the expectation value of n given in Eq. (13T1) of the Appendix 
we obtain the energy surface 



^(0 



A+jA 
2j 



A 



27 



c 



L 



W4^i(-c^)A2r(-c^) 



-J + 1 < A < J 



(23) 



2i ^2i--lV S J/ 

and 7^(0 = -|(1 - A) for A = -j. 

We calculate this expression as a function of A and 7, and for all possible values of A with 
A G [—2, 4] and 7 G [—5, 5], which allows us to choose the value of A for which the energy is 
minimum and onto which the coherent state is to be projected. By noting that the classical 
value Ac of the constant of motion, given in Ec. ([6]), when rounded to the nearest integer or 
half-integer, never differs by more than one unit from its exact quantum counterpart, the 
minimizing procedure simplifies to testing for |A — Ad < 1, and approximating A to the 
nearest integer (if is even) or half- integer (if is odd). The minimum energy so obtained 
is plotted in Fig. ([2]), while the resulting constant of motion and its contour levels are shown 
in Fig.lE]). 

While similar, the behaviour for A > 1 appears to be different than that for A < 1. 
The energy surface shows a fictitious phase transition near 7 ^ 0.7 for A > 1 (cf. right of 
Fig. ([2])) that does not exist for A < 1. In fact, the regime A > 1 appears only when < 0; 
one may visualize physical situations where this may occur: if the atoms are immersed in an 
external magnetic field and the energy levels appear as a Zeeman splitting of spectral lines, 
one may think of continuously tuning the field intensity until a sign reversal is obtained, thus 
interchanging the excited and base levels. It will be seen, however, that the exact quantum 
solution does not show this transition (cf. Fig. ([6]) below). The case A = 1 (or equivalently 
uja = 0) is, nevertheless, a special and boundary case: the Hamiltonian in Ec. ([T]) simplifies 
considerably, the atomic levels are degenerate, the quantum energy spectrum (see Section 
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FIG. 2: Energy surface of the ground state used to calculate the value A of the constant of motion, 
for = 20 atoms. The zoom at right shows phase transitions when A is greater than 1 (see text 
for details). 

V) is also highly degenerate, and the poles (both North and South) contract to a single 
point. In fact (cf. Ec. ([6])), the North Pole is no longer a minimum critical region, but a 
saddle point. 
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FIG. 3: Constant of motion Xmin determined by minimizing the expectation value of the Hamil- 
tonian with respect to the projected state, as a function of the interaction strength 7 and the 
detuning parameter A, for = 20 atoms. Contour levels are shown at right. 

From the expressions derived for the observables in Section III, their expectation values 
can also be plotted as functions of the interaction strength 7 and the detuning parameter A. 
It will be seen that their behavior will be inherited from that of Xmin itself. As an example, 
we show {Jz) in Fig.(jll). 

There is no field squeezing since, from the quadrature components, (g) = (p) = and 
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FIG. 4: Expectation value of the population inversion operator J^, for the variational ground state, 
as a function of the interaction strength 7 and the detuning parameter A, for = 20 atoms. 



(g^) = (n) + h = {fP'). However, for the matter squeezing coefficient we have [32| 



'2(AJi 



(24) 



where J± is a component of J transverse to (J); as (J) = {Jz)e.z, we take Jx to be this 
component, and show ^ in Fig. ([5]) (left). ^ gives us information of how good the trial 
projected state is approximating the exact solution of the eigenvalue problem of the TCM 
Hamiltonian, and this comparison will be made in Section V below. 





FIG. 5: Squeezing parameter ^ for the matter (left), and entanglement entropy Se between field 
and matter (right), for the variational ground state, as a function of the interaction strength 7 and 
the detuning parameter A, for = 20 atoms. 



Finally, the entanglement entropy 5*^; is zero for coherent states since these are expressed 
as product states. In [26|, in order to calculate a useful expression for Se, we first traced 
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over the field and then again over one of the matter modes. When a projected state is used, 
it is no longer a product state and we may trace over the matter (or equivalently over the 
field) to obtain 

^i? = -$^7^.1ogn (25) 

V 

where Vv is already normalized, as given in Ec. (fT4l) . This is shown in Fig. ([5]) (right). Its 
behavior is very close to that found in [26] . 



V. PROJECTED COHERENT STATE VS. EXACT SOLUTION 



The exact solution for the ground state in the Tavis-Cummings model was presented 



m 



261], where the natural basis \v) ® |j, \ — v) was used (natural because the Hamiltonian 



has A as a constant of motion): 

A+i 

= 5^ \v) ® |j, \ -v) , (26) 

j^=max[0,A— j] 

Analytical solutions for values of A up to A = — y + 4 were also given, both in resonance 
(A = 0) [33I] and away from resonance (A 7^ 0) 26]. For greater values, numerical solutions 
were obtained. 

When the exact solution is compared with that obtained from the projected state, we 
find that the quantum phase transitions are exactly reproduced, as shown in Fig.(l6l), even 
for a large number of atoms (A^ = 20) and when away from resonance. The straight lines 
in the lower left figure correspond to the energy of the ground state for different values of 
the constant of motion A, starting from —10 (0 photons) and up to —6 (4 photons). In each 
case the eigenstate is a mixture of to A + ^ photons. The projected state approximation is 
stunningly close; so much so, that both graphs (quantum and projected state) practically lie 
on top of each other and are indistinguishable. Whereas the projected state has a constant 
value of the ground state energy inside the poles (North (|0) ® |j, — j)), and South (|0) ® 
|j, j))), the exact quantum solution does not when A > 1 {uja < 0); the solution therefore 
is very precise for a detuning of A = 0.2, as shown in the left side of the figure, but fails 
slightly near the boundary of the poles for A = 1.5, as shown at lower right. 

The case A = 1 is a particular boundary, as mentioned in the previous section, since it 
simplifies the Hamiltonian ([3]) making it independent of the inversion population operator 
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FIG. 6: Ground state energy per particle of the projected state solution compared with the exact 
quantum solution, as a function of the interaction strength 7 for = 20 atoms. The plots on 
the left correspond to a detuning parameter A = 0.2, while those on the right to a A = 1.5 (and 
therefore loa < 0). The bottom plots are close-ups of those on top. Note that both solutions lie on 
top of each other and are indistinguishable, except for the close-up for A = 1.5. 

Jz- In this case, if 7 = also, we have a dense degeneracy of the Hamiltonian eigenvalues: 
the matter does not see the radiation field, so the degeneracy is that of the atoms themselves, 
i.e., 2j + 1. Any small deviation of 7 away from zero unfolds the energy levels, as shown 
in Fig. ([7]). In this case our method for choosing the A- value onto which to project fails, as 
2j + 1 values yield the same energy. We shall, in what follows, choose values for A above 
and below 1 to exemplify the behavior of observables. 

The constant of motion, Xmin is shown in Fig. (JHl) as a function of the interaction parameter 
7 for both, the semiclassical projected and the quantum cases. Its discrete behavior is 
inherited by all the quantum observables of interest. The slope observed in the A-steps is 
given precisely by the difference in energy between the absorbed photons and the atom's 
energy level separation. This is a consequence of the system being away from resonance 
(A 7^ 0). When in resonance, the steps are horizontal. The plots at left show A for a 
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FIG. 7: Energy spectrum for A = 1. When 7 = we have a (2j -1- l)-fold degeneracy (left plot 
shows the first 2 such levels), which unfolds as soon as 7 is perturbed away from zero. Spectra for 
7 = 0, 0.1, 0.2, and 0.3 are shown. 

detuning of A = 0.2, while those at right show it for A = 1.5. Note (lower right plot) that 
even when Ua < the projected state solution tries to reproduce the non-constant behavior 
near the pole boundary. Otherwise the graphs are indistinguishable. 
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FIG. 8: Constant of motion Xmin of the projected state solution compared with the exact quantum 
solution, as a function of the interaction strength 7 for = 20 atoms. Note, once again, that 
both graphs practically lie on top of each other and are indistinguishable. The plots on the left 
correspond to a detuning parameter A = 0.2, while those on the right to a A = 1.5 (and therefore 
UJA < 0). The bottom plots are close-ups of those on top. 



The photon number fluctuations (An)^, as noted in the Introduction, are very well re- 
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sembled by the projected state solution to the TCM (cf. right of Fig.([T])). The differences 
are minute. The same trend is found for the dispersion in J^, as is to be expected, since 
(J^) = A - (n), (J2) = (n^) - 2\{n) + A^, and therefore (AJ^)^ = (An)^. Both, the expec- 
tation value for and its dispersion are plotted in Fig. ([2]) as functions of the interaction 
strength 7, slightly away from resonance (A = 0.2). Note the scale in the ordinate axis for 
(AJ,)2. 
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FIG. 9: Expectation value for and its dispersion (AJ^)^ of the projected state solution compared 
with the exact quantum solution, as functions of the interaction strength 7 for = 20 atoms. Both 
graphs practically lie on top of each other (cf. the scale for (AJ^)^). The plots correspond to a 
detuning parameter A = 0.2. 



As a signature of the goodness of our trial projected state in reproducing the exact 
quantum solution, one may use the behavior of the squeezing spin coefficient ^ as given by 
Eq. ([23]). Even though it greatly exceeds the value of 1 away from the poles (cf. Fig. flTUl) ). 
thus suggesting a small overlap with the quantum state, both results lie exactly on top of 
each other. Once again, for a large detuning making tu^ < 0, the round shape of the pole is 
badly approximated at its boundary. For A < 1, however, the approximation is exact. 

By taking the trace with respect to the field (matter) states the reduced density matrix 
takes the form 



^matter 



field 



mm{A+-f,7V} 

E 

n=0 
A+j 

;/=max{0,A— j} 



— n){N — n, n\ 



(27) 
(28) 
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FIG. 10: Squeezing parameter ^ of the projected state solution compared with the exact quantum 
solution, as functions of the interaction strength 7 for = 20 atoms. The plots correspond to 
a detuning parameter A = 0.2 (left) and A = 1.5 (right). The round shape of the pole is badly 
approximated at its boundary, when A > 1, while the approximation is exact for A < 1. 



where c^j^n_^ (or c,y) is determined from the Hamiltonian diagonalization. As the reduced 
density matrix of the matter is diagonal, the matter-field entanglement entropy equals that 
between the atoms occupying the two hyperfine levels 



mm{X+f,N} 



- E 



n=0 



In 



(29) 



in both cases. This may be compared with that for the projected state, Ec. (!25|) . to give 
the result shown in Fig. ffTT]) . The same behavior of excellent approximation for A < 1, and 
poor approximation at the south pole boundary when A > 1, is obtained. 



VI. DISCUSSION AND CONCLUSIONS 

Our approximation is not exact: several important observables such as Se, (Ag)^, {n), 
and (An)^ are symmetric under a reflection about A = 1 for the projected state, while this 
cannot be true for any observable coming from the Hamiltonian ([3]); as an illustration of 
this we show (An)^/A^^ as a function of A, for both the projected state (light curve) and 
the quantum state (dark curve), in Fig. flT^ . for 7 = 0.75. While the projected state has a 
constant value of at both poles, the quantum state decays to asymptotically as it enters 
the South Pole. 

A good measure of the distance between quantum mechanical states is given by the 
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FIG. 11: Entanglement entropy Se of the projected state solution compared with the exact quan- 
tum solution, as functions of the interaction strength 7 for = 20 atoms. The plots correspond 
to a detuning parameter A = 0.2 (left) and A = 1.5 (right). The round shape of the pole is badly 
approximated at its boundary, when A > 1, while the approximation is excellent for A < 1. 
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FIG. 12: (An)^ seen as a function of of A for the classical projected state (light curve), compared 
with the exact quantum solution (dark curve), for = 20 atoms and 7 = 0.75. 

fidelity; for pure quantum states it measures their distinguishability in the sense of statistical 
distance 27|, but it is customary to use the fidelity as a transition probability regardless of 
whether the states are pure or not. Measuring the fidelity between our projected state and 
the exact quantum ground states 



(30) 



gives a result very close to 1, except in the region of classical phase transitions. Fig. (fT3|) 
(left) shows the result as a function of 7 for a detuning parameter of A = 0.2. F = 1 
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inside the North Pole and drops to F = 0.996 when crossing the separatrix into the Parallels 
region, only to approach its value of 1 again as I7I continues to increase. We have seen that 
the South Pole is less well represented by the projected state for A > 1; the behavior of 
F as a function of A near the South Pole is shown in Fig. flTSl) (right), for 7 = 0.75, for 
which value the South Pole lies at A > 1.56. It is seen that even in this case F drops to 
approximately 0.6 and quickly recovers. 




FIG. 13: Fidelity F between projected state and exact quantum ground states, as a function of 
the interaction strength 7 (left) and of the detuning parameter A (right), for = 20 atoms. The 
graph on the left corresponds to a detuning parameter A = 0.2, and shows F = 1 inside the North 
Pole dropping to F = 0.996 when crossing the separatrix into the Parallels region. The one on the 
right corresponds to an interaction strength 7 = 0.75, showing F ~ 0.6 at the South Pole. 

It has been emphasized that the Dicke model presents quantum phase transitions when 
ujp'-f = Though difficult to satisfy for optical systems, some proposals to overcome 

the experimental problems have been reported |34|. 

More recently jssi, the Tavis-Cummings Hamiltonian has been physically realized for 
artificial atoms in the form of superconducting q-bits at fixed positions and coupled to a 
resonant cavity mode. The predicted -\/iV-behavior for the collective A^-atom interaction 
strength of the model is observed experimentally in good agreement. It is argued that the 
presented approach may enable novel investigations of superradiant and subradiant states 
of articial atoms. 

In this work we have shown that an excellent approximation to the exact quantum solution 
of the ground state of the Tavis-Cummings model is obtained by means of a semi-classical 
projected state. This state has an analytical form in terms of the model parameters and 
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allows for the analytical calculation of the expectation values of field and matter observables, 
entanglement entropy between field and matter, and squeezing parameter. In our discussions 
we have taken the picture of two-level atoms interacting with a single mode electromagnetic 
field, but by re-interpreting the parameters the results are useful to describe, for example, 
the superradiant properties of condensed matter phases, or linear ion traps proposed for 
quantum computation. 

Acknowledgments 

This work was partially supported by CONACyT-Mexico and DGAPA-UNAM. 



Appendix A. Expectation Values of Field Operators 



The expectation value of the photon number operator can be found analytically using 
the overlap f|T3l) . with (' = (, to get 



in) 



iX + j)-i^+j)L',-^ii-v)/L',-\-v) , A>j 



(31) 



In the same way we find that the squared fluctuations are 



(Ah) 



j- 



-XI 



-2j 



-XI 



(A + j) {L',;I,{-x)/l',;\-x) + (A + J - 1) L',T^,{-x) / L',-\-x) 



-(A+J) 



L',;l,{-x) L',-\-x) 



-J + 1 < A < i 



A>J 
(32) 
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